Power-law distributions of particle concentration in free-surface flows 
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Particles floating on the surface of a turbulent incompressible fluid accumulate along string-like 
structures, while leaving large regions of the flow domain empty. This is reflected experimentally by 
a very peaked probability distribution function of c r , the coarse-grained particle concentration at 
scale r, around c r = 0, with a power-law decay over two decades of c r , If(c r ) oc c~ l3r . The positive 
exponent /3 r decreases with scale in the inertial range, and stays approximately constant in the 
dissipative range, thus indicating a qualitative difference between the dissipative and the inertial 
ranges of scales, also visible in the first moment of c r . 

PACS numbers: 47.27.-i, 47.27.od, 47.52. +j 
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1. INTRODUCTION 

Consider the introduction of passive, neutrally buoy- 
ant tracer particles into an incompressible turbulent flow, 
for example ink injected into a vigorously stirred tank of 
water. After a transient mixing period, the dye becomes 
uniformly distributed. Further stirring will have no ob- 
vious effect; the dye particles will be rearranged on a 
microscopic scale, but on larger scales, the concentration 
distribution will remain uniform. 

Now consider a different experiment in which the con- 
taminant particles are lighter than the host fluid, so that 
they float on the fluid-air interface. This situation, which 
occurs naturally for certain classes of phytoplankton or 
contaminants [ljl, 0], has been studied only recently in 
the laboratory 4[. Even when the floaters are initially 
distributed uniformly on the surface of the host fluid, 
which is continuously stirred before and after the surface 
is covered, the concentration distribution of the floaters 
will not remain uniform. Rather, the floaters will clus- 
ter into string-like structures, see FigfTJ The underlying 
turbulence, which drives the motion of the floaters, will 
cause the string-like distribution to fluctuate in space and 
time, but there will be no tendency toward concentration 
uniformity. 

This clustering is not caused by an interaction between 
the floaters, nor is it produced by surface tension effects 
or by wave motion, ft is seen in simulations of the tur- 
bulent motion of the underlying fluid when these compli- 
cating effects are not included. In the simulations 
the fluid, assumed to be governed by the incompressible 
Navier-Stokes equation, is driven at large spatial scales 
to turbulence. The floaters merely follow the fluid mo- 
tion and sample the horizontal velocity components at 
the surface; they are therefore treated as inertia-free. 

Here we study the Lagrangian evolution of the steady- 
state tracer concentration of the floaters at the interface 
between air and water. Although the underlying turbu- 



lence is incompressible, the particles which are lighter 
than water are constrained to move only along the two- 
dimensional interface, they cannot follow water into the 
bulk. The tracers flee fluid up-wellings (sources) and 
cluster along the fluid down-wellings (sinks) and form 
a compressible system. From a dynamical systems per- 
spective the effect seen in Fig. [T] results from a com- 
petition between stretching along one spatial direction, 
and contraction over another one. The particle motion 
is characterized by a positive and a negative Lyapunov 
exponents, with the sum being negative @, 0|- 

The effect of the surface flow that induces floater clus- 
tering is quantified by the two-dimensional compress- 
ibility: C = (|(V2 ■ V ^ \ / (ty?,v^ V For isotropic 

turbulence, C lies between (incompressible fluids) and 
unity (potential flows). Past experiments and Direct Nu- 
merical Simulations (DNS) have measured C w 1/2 0], 
making this a strongly compressible system. The accu- 
mulation of particles in surface flows is reminiscent of the 
physics leading to clustering of inertial particles in incom- 
pressible flows [H, [H . There the clustering arises due to 
inertia of advected particles whose trajectories deviate 
away from the stream lines. In the present case of free- 
surface turbulence, the floaters are passive and closely 
follow the local flow field, the clustering arises purely 
because their motion is constrained to the air-water in- 
terface. The fluid flow acts as a dynamical system on the 
floating particles [!, @, H, [l(J U , 13] ■ At scales below the 
dissipative scale of the turbulence, rj, the flow is viscous 
(smooth), and theories based on the stretching of fluid 
elements lead to the prediction that the particle distri- 
bution is multifractal 0, [T3| ■ For particle distributions in 
the inertial range where the velocity field is not smooth, 
the only available study comes from a surface flow model 
that suggests a multi-fractal particle distribution for both 
the inertial and dissipative ranges Q. 

We study the concentration statistics of particle clus- 
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ters in free-surface flows as a function of the dimension- 
less coarse-graining scale r in the steady state. The dis- 
tribution of the coarse-grained concentration n r over a 
domain of scale r, is very intermittent (Fig. [I}, the more 
so as the scale r decreases. This can be seen in Fig. [TJ 
where particles have been expelled from large regions of 
the flow by surface sources and have accumulated along 
dense line-like surface sinks. Intermittency can be quan- 
tified by measuring the moments of n r averaged over all 
points in space (Eulerian representation) . In the present 
work, we instead investigate the Lagrangian concentra- 
tion c r averaged over a domain of size ~ r around each 
particle, which is equivalent to the Eulerian quantity (see 
Eq. H}. 

The first main result of this work is the observation of 
two different scale- free regions for the first moment of c r , 
corresponding to the dissipative and inertial ranges, in 
qualitative agreement with 0]. The PDF of the coarse- 
grained concentration Tl(c r ) exhibits a steep (discussed 
below) power-law behavior at small values of c r , char- 
acterizing the regions with a very low particle density. 
PDFs with power-law behavior are frequently encoun- 
tered in nature in many different contexts [13| . usually 
over a limited range of values of the random variable. We 
find that the PDF evolves as a function of scale r, and 
demonstrate that this evolution is qualitatively different 
in the dissipative and inertial ranges. 



2. EXPERIMENT 

The experiments were carried out in a tank lm x lm 
in lateral dimensions. It was filled with water to a depth 
of 30cm. Turbulence was generated by a pump (8hp) 
which re-circulates water through a system of 36 rotating 
jets placed horizontally across the tank floor. This sys- 
tem ensured that the source of turbulent injection was 
far removed from the free-surface where the measure- 
ments were made More importantly, it minimized 
the amplitude of surface waves which were unavoidable, 
but have been proven to play no role in the effects re- 
ported here 0]. The hydrophillic particles chosen here 
are subject to capillary forces which are very small com- 
pared to forces coming from the turbulence, and do not 
affect the results as they do in [ToL 14 1. The non-inertial 



character of the particles is minimal because the product 
t s Ai, where Ai is the largest Lyapunov exponent and t s 
is the stopping time of the particle [l2j , is of the order 
5 x 10 -4 , which is too small to lead to any significant 
inertial effect 
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During an experimental run, particles (50fj,m diameter 
and specific gravity of 0.25) were constantly seeded into 
the fluid from the tank floor. From there they under- 
went turbulent mixing as they rose due to buoyancy and 
were uniformly dispersed by the time they rose to the 
surface. Once at the free-surface, their motion was con- 
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FIG. 1: (Color online) Top panel: Visualization of the steady- 
state particle distribution shows string-like coagulations sur- 
rounded by (relatively) depleted regions. Bottom panel: The 
Eulerian concentration field n r in the steady-state shows a 
very intermittent distribution, with very large spikes concen- 
trated along string-like structures, surrounded by vast de- 
pleted regions which are nearly flat. The x and y axis are 
in units of r\ (see text). The coarse- graining scale chosen here 
is r = 3.4. Decreasing the coarse-graining scale r increases 
the intermittent nature of the particle concentration. 



strained to the two-dimensional surface plane and they 
were tracked via a high-speed camera (Phantom v. 5) sit- 
uated above the tank. Instantaneous velocity fields were 
measured using an in-house developed particle imaging 
velocimetry (PIV) program which processed the recorded 
images of the floaters. The scheme of particle injection 
ensured that both the sources and sinks of the flow re- 
ceived equal coverage of particles to perform the PIV 
analysis. The injection scheme was also necessary to re- 
place particles that left the camera field of view, allowing 
the experiment to be performed for an adequate time pe- 
riod. The experimentally measured velocity vectors were 
spaced (on average) by Sx = 2.5 r\. The level of tur- 
bulence achieved with our setup leads to a Taylor Mi- 
croscale Reynolds number Re\ in the range 150 — 170. 
In the flow, the integral length scale is Zo ~ 1.45cm, the 
Taylor scale is A f=a 0.45cm and the dissipative scale is 
V — (^ 3 /e) 1 { 4 ~ 0.02cm, where e is the energy dissipa- 
tion rate [10]. The camera field-of-view is a square area 
of side length L = 9cm. The camera's height above the 
water surface was chosen so that a pixel size is roughly 
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0.1mm, comparable to the dissipative scale of the turbu- 
lence. 

The measured velocity field was then used to solve 
the equation of motion for Lagrangian particles : = 
v(xj(t),t), where v(xj,i) is the velocity field and Xj = 
( x iTUi) are the individual particle positions. To achieve 
sub-grid resolution of the measured velocity field, as is re- 
quired to measure n r in the dissipative range, the bi-cubic 
interpolation scheme discussed in [13] was implemented. 
To use this scheme it is necessary for the measured veloc- 
ity grid spacing to satisfy the criterion 5x < ttt], where 
Sx is the (above quoted) average measured velocity grid 
spacing. We have tested to ensure that the results pre- 
sented here do not depend on the velocity grid spacing 
by varying the spacing from 5x = 2.5 77 to 5x = 4 77. 
All of the statistics presented in this work were obtained 
by evolving a total number of N t ~ 4 x 10 5 Lagrangian 
particles per frame for each experimental run. For each 
data set, a uniform distribution of Lagrangian tracers was 
generated at time t = 0s and evolved using the above 
described scheme. The statistics presented here are dis- 
cussed in terms of length s compared to the dissipative 
length scale 77, r = s/rj, making r dimcnsionlcss. 

Concentration statistics in turbulent studies can be 
studied by dividing the field of view into boxes of size 
r and counting in each cell the number JV(r, t) of parti- 
cles, thus computing the the value of the corresponding 
Eulerian (coarse-grained) concentration, n r (t) = (^^t)) ■ 
Alternatively, the concentration can be measured in the 
Lagrangian frame. To do this, one draws a box (or cir- 
cle) , of size r around each (moving) particle and counts 
the number of particles N*(r,i) that lie inside, including 
the centered particle. The Lagrangian (dimcnsionlcss) 
concentration is then defined as: c r (t) = ^ N ^t)) ' wnere 
the denominator is the Eulerian frame mean. Physically, 
averaging a quantity over a set of Lagrangian particles is 
equivalent to averaging over space, with a weight propor- 
tional to the particle density. The relation between the 
moments of n r and c r : 



(c P r) = W +1 >- 



(1) 



follows from this simple fact, at least when the value of r 
is small enough, so the correlation functions of the con- 
centration field simply scale with r (is| . Likewise the 
PDF's of n r and c r are related: ra r II(ra r ) = n(c r ). Ex- 
perimentally, one obtains much more accurate results by 
using the Lagrangian field c r , rather than n r filj ]. We 
have checked that the accurate results determined with 
c r are consistent with the results obtained with n r , using 
Eq[T] In the following, we express our results in terms 
of the Lagrangian concentration field, cy. The statistics 
of the concentration field discussed here have been taken 
only in the statistically steady-state regime, by systemat- 
ically leaving out the initial transients, necessary for the 



TABLE I: Turbulent parameters measured at the surface. 
Measurements are made at several values of the Re\ with 
an average Re\ ~ 160. The parameters listed are averages, 
with deviations less than 10%. 



Taylor microscale 


V ((s^/e-) 2 > 


U.4i (cm) 


Taylor Re a 


Re v rms X 

V 


160 


Integral Scale 


, n , (v\\ (x-\-r)v\\ (x)) 


1.42 (cm) 


Large Eddy Turnover Time 


^ v rm.p 


0.43 (s) 


Dissipation Rate 




6.05 

(cm2/s3) 


Kolmogorov scale 




0.02 (cm) 


RMS Velocity 


Vrms — J(V 2 ) - (v) 


3.3 (cm/s) 


Compressibility 




0.49 ± 0.02 



moments of the distributions to reach their steady-state. 



3. RESULTS 

We discuss the properties of the steady-state distribu- 
tion of the coarse-grained particle concentration c r de- 
fined above. We begin with the lowest non trivial mo- 
ment : (c r ), and its dependence on the dimcnsionlcss 



coarse-graining scale r [12|]. Scaling relations, such as 
(c r ) ~ r~ a , provides a characterization of the topologi- 
cal structure of the particle distribution [l9[ ■ For simple 
geometrical structures such as points, lines, and surfaces, 
the scaling exponent a take on the simple integer values 
a = 2,1,0, respectively. For more complicated struc- 
tures, as in our case, a is a non-integer and gives a quanti- 
tative measure of the complex topological structure. Fig. 
[5] shows that (c r ) exhibits two different power-law be- 
haviors, at small (r < 3) and at larger (r > 5) scales. 
These data can be fit with an exponent = 0.92 ± 0.02 
{pci = 0.79 ± 0.03) in the inertial (dissipative) ranges. 
According to Fig. [21 the dissipation range extends up 
to scales larger than 77 (for r « 3 — 5). This is consis- 
tent with the fact that the dissipative range crossover in 
many experiments occurs at r > 1 [20I [21| . The differ- 
ence between the two scaling domains can be best seen 
by plotting the derivative dln(c r )/dlnr, see the inset in 
Fig. The qualitative impression, see FigO] that the net 
result of the coagulation process are line-like structures 
is completely consistent with the values of the exponents 
close to 1. Higher order moments of the coarse-grained 
concentration also exhibit unique scaling regimes in the 
inertial and d issipative ranges, transitioning between the 
two regimes at approximately brj. Our results are thus 
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FIG. 2: (Color online) First moment of the Lagrangian con- 
centration field c r . The average (c r ) exhibits two scaling 
regimes. In the dissipative range, for s < 3n, (c r ) scales as 
r ~ a d w ith oid = 0.92 ± 0.02, whereas for s > 5n, it scales as 
r~~ a ' with on = 0.79 ± 0.03. The inset shows the derivative 
dln d ^ T ' > , and demonstrates the quality of the reported power- 
law behavior. 

in qualitative agreement with the observations of 0] for 
(njr), which also show two separate scaling intervals cor- 
responding to the inertial and dissipative ranges. 

We now discuss the PDF II(cy) of the concentration. 
Fig. [3] shows the PDFs of c r , normalized by the mean 
(c r ), in the dissipative and inertial ranges. In the iner- 
tial range, and for r = 3.4 in the dissipative range, the 
PDFs display a clear maximum in the distribution, with 
this maximum occuring at smaller values of c r j (cy) as 
r decreases. For values of c r /{c r ) less than this max- 
imum, the PDFs exhibit a plausibly linear behavior : 
II(c r ) oc c r , which corresponds to a PDF in the Eulerian 
frame H(n r ) oc constant (see Eq. ([1])). This ensures that 
the PDF is normalizable in both the Lagrangian and Eu- 
lerian frames. This very small c r /(c r ) wing of the PDF 
is not observed deep in the dissipative range, for r < 1.1, 
due to the fact that values of c r /{c r ) corresponding to 
less than one particle are not adequately resolved. 

The most interesting aspect of the PDFs in the dissi- 
pative (inertial) range is the power-law behavior at low 
(intermediate) values of cy/ (c r ) greater than the PDF's 
maximum. This algebraic decay corresponds to the de- 
pleted regions seen in Fig. Q] rather than the string-like 
coagulations. In both spatial ranges, the PDFs can be 
characterized by a power-law behavior, II(c, ) oc c~^ r , 
over a limited range of c r /(c r ) followed by a faster than 
algebraic fall-off at large c r /(c r ). The variation of the 
exponent f3 r in the dissipative range, Fig[3l is very small 
: fl r ~ 0.8 ± 0.05 for r < 3.4. In comparison, the ex- 
ponent (3 r decreases when the scale r increases in the 
inertial range (see Fig. [3]), with (3 r ~ 0.5 at r = 28.6 and 
(3 r ~ 0.75 at r = 5.0. The fast decay of the PDFs at large 
c r /{c r ) ensures that all the moments of the distribution 
(c£) exist, for all p > 1. 

The PDFs in both the inertial and dissipative ranges 
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FIG. 3: (Color online) Top panel : II(c r ) in the dissipative 
range r < 3. The PDFs of the concentration, scaled by its 
mean value (c r ) exhibit a power-law distribution at low value 
of c r /(c r ): R(x) oc x~> 3 with /3 ~ 0.8 ± 0.05. The PDFs 
decay faster than algebraically for c r /(c r ) > 1. Bottom panel 
: II(c r ) in the inertial range r > 5. The PDFs in the inertial 
range also exhibit a power-law distribution at low value of 
c r /(c r ) greater than the PDF's maximum : n(c r ) oc c^T 
whose exponent /3 r decreases when r increases. Here, f5 r ~ 0.5 
at r — 28.6 and increases to f) r ~ 0.75 at r — 5.0 The decay of 
the PDF at large values of c, /(c r ) > 1 is faster than algebraic, 
clearly seen for r = 7.9 and 5.0. The PDFs in the inertial 
range have been shifted vertically to better see the evolution 
of the power-law region. 

show a systematic variation of their form as a function of 
scale r, qualitatively similar to the evolution of the veloc- 
ity increment PDF in turbulence (ljj . In the dissipative 
range, this observation is consistent with the theoreti- 
cal arguments of @, leading to the prediction that the 
particle distribution is multi-fractal. Such a multi-fractal 
property implies that the PDF must depend on the scale 
r. The change in the PDF with scale r mostly affects the 
large c r j (c r ) behavior, where the PDF decays faster than 
algebraically. The pre-factor of the algebraic behavior of 
the PDF also varies with scale r. 

4. SUMMARY 

To summarize, we have characterized the statistical 
properties of particle clustering on the surface of a turbu- 
lent flow by measuring the coarse-grained concentration 
around each particle in the system. Our study allows us 
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to cover a range of scales extending from dissipative up 
to inertial scales. The first moment of the particle distri- 
bution (c r ) exhibits two different scaling regimes in the 
dissipative and in the inertial ranges, for scales smaller 
and larger than ~ 5r/. The PDF of c r vary systemati- 
cally with scale. Indeed, our limited results are consis- 
tent with a multi-fractal scaling both in the dissipative 
and in the inertial ranges, which is generally consistent 
with the findings of [3] • The PDFs exhibit power-law be- 
havior at small values of c r /(c r ), characterizing the very 
high probability of having regions containing very low 
concentrations. 

The corresponding distribution of the Eulerian con- 
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centration n r behaves as II(n r ) oc n r 



(l+/3r) 



which is not 



even normalizable when n r — > 0, thus implying satura- 
tion at a small value of Tip . £lS indeed su gges ted by Fig. [3] 
Such a behavior has been observed in 11 01 with /3 r = 1. 
The explanation in terms of caustics [22j proposed in 
is unlikely to apply here, both because our particles 
strictly follow the Lagrangian evolution equation, and 
because (3 r ^ 1. 

Our work thus indicates that the intermittent distribu- 
tion of particles in free-surface flows is characterized both 
by very high particle concentration along string-like re- 
gions as well as by a strong depletion of particles in other 
parts. A proper description of the resulting power-law 
tails requires a better understanding of the very efficient 
expulsion mechanism of particles from large regions of 
the flow. Understanding this particle clustering is of fun- 
damental interest and practical use in a range of surface 
transport phenomena, such as the dispersion of phyto- 
plankton or contaminants in oceanic flow P, 0, 01 ■ 
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